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Abstract. Transport in suspended metallic single wall carbon nanotubes in the 
presence of strong electron-electron interaction is investigated. We consider a tube of 
finite length and discuss the effects of the coupling of the electrons to the deformation 
potential associated to the acoustic stretching and breathing modes. Treating the 
interacting electrons within the framework of the Luttinger liquid model, the low- 
energy spectrum of the coupled electron-phonon system is evaluated. The discreteness 
of the spectrum is reflected in the differential conductance which, as a function of 
the applied bias voltage, exhibits three distinct families of peaks. The height of the 
phonon-assisted peaks is very sensitive to the parameters. The phonon peaks are best 
observed when the system is close to the Wentzel-Bardeen singularity. 
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1. Introduction 

Carbon nanotubes are new prototypical candidates for one-dimensional nanoconductors. 
As such, electron transport through such macromolecules has been investigated 
extensively. Single wall carbon nanotubes (SWNTs) behave as metallic (n — m = 31, I 
is an integer) or semiconducting conductors (n — m ^ 31) depending on their wrapping 
vector, (n,m) = na[ + ma^, where ai, 02 are the graphene primitive lattice vectors. 
PUI21E]- Moreover, Coulomb blockade phenomena, the Kondo effect, and Tomonaga- 
Luttinger liquid properties have been experimentally demonstrated [U 03 El El IE] • Such 
measurements confirm the relevance of electron correlations effects in SWNTs at low 
energies as predicted in earlier theoretical studies HH1 HI]- Recently, suspended 
nanotubes have been fabricated, and the transport properties have been measured by 
several groups [121 EH EH EH EE]- Since the suspended nanotube is free from the 
substrate, it is expected that it can easily vibrate [THl HZI- At the same time, some 
theoretical works have started to investigate the effects of the mechanical motion on 
the electronic transport properties of interacting SWNTs [TBI EH 1201 121] • While Ref . 
[T9l l2T] focus on Coulomb blockade phenomena, the SWNT is treated as a Tomonaga- 
Luttinger liquid |22] with periodic boundary conditions (valid for very long nanotubes) 
in Ref. JHH201- However, typical SWNT lengths are of the order of sub-micron to micro- 
meters, yielding a mean level spacing for the electronic energy levels in the meV energy 
range. This discreteness of the energy levels can be clearly seen in low temperature 
measurements, e.g. by measuring the excitation spectrum of a SWNT quantum dot 
[2*3*] . In finite length SWNTs, the quantized nature of the vibrational modes of a SWNT 
should be considered as well at low enough energies. 

In this paper, we investigate the effects of the mechanical vibrations on the 
electronic properties of a doubly clamped SWNT of finite length within a full quantum 
mechanical treatment of the low energy excitations. Even if the metallic condition 
n — m = 31 is satisfied, small energy gaps ~ lOmeV appear due to curvature effects. The 
armchair nanotubes (n = m), however, still stay metallic because of their high symmetry 
[21] . Therefore we consider the armchair SWNTs in this paper. The interacting electrons 
in the nanotubes are considered within the standard bosonization framework of the 
Tomonaga-Luttinger liquid model with open boundary conditions [213 ES]- where the 
electronic excitations (charged and neutral plasmons) have bosonic character. The 
low energy vibrations of the nanotube can be described in terms of acoustic phonons, 
associated to the longitudinal stretching mode and to the transverse twisting and 
breathing modes [22j. The coupling between electron and phonons is described in terms 
of the deformation potential caused by the acoustic modes. In our model this yields 
a bilinear bosonic coupling between phonons and plasmons. In Ref. [20] , the acoustic 
modes have been treated within a continuous elastic model, and have been integrated 
out in order to obtain a phonon-mediated retarded electron-electron interaction. It was 
shown that while the twisting mode can be neglected at low energies, the effect of the 
breathing mode and stretching modes should be considered. In particular, retardation 
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effects of the breathing mode can be neglected, and an effective renormalization of 
the electron-electron interaction strength is found with possibly attractive interactions. 
In contrast, retardation effects due to the stretching mode must be retained, so that, 
due to the combined effects of stretching and breathing modes, the Wentzel-Bardeen 
singularity, describing a superconducting instability of the electron-phonon system [28J, 
can be reached. In this paper, we consider the case in which electron as well as phonon 
modes are quantized, and look to the energy spectrum of the interacting electron-phonon 
system. As in Ref. [20], the major effect of the breathing mode is to renormalize the 
electron-electron interaction strength. On the other hand the electronic as well as the 
stretching-phonons parameters get renormalized by the mutual interaction, yielding 
"electron" -like and "stretching-phonons" -like excitations. Such electron-like excited 
levels appear as large peaks in the differential conductance, while stretching-phonon- 
like assisted processes appear as small peaks beside the large peaks. The phonon peak 
height is very sensitive to the system parameters, and it becomes larger and larger when 
the value of the charged plasmons' velocity is close to the sound velocity, i.e. when the 
Wentzel-Bardeen instability [2E| is approached. Thus, the experimental observation of 
the phonon peak would be an indication of the strong fluctuation of the superconducting 
order parameter, which corresponds to the Wentzel-Bardeen singularity, in SWNTs. At 
temperatures higher than the level spacing, the discreteness of the peaks is washed out. 
The differential conductance then shows the characteristic Tomonaga-Luttinger power- 
law dependence on voltage, with exponents modified by the electron-phonon interaction. 

2. Low-energy Hamiltonian 

To start with, we consider a suspended, doubly clamped metallic SWNT whose fixed 
ends are at x = and x = L. We wish to develop a low energy theory for suspended 
metallic SWNTs, and restrict ourselves to the two lowest linear in k electronic subbands. 
Moreover, in experiments the nanotubes are typically intrinsically £>-doped such that the 
Fermi level is shifted away from half-filling. Long-range Coulomb interaction is dominant 
in isolated SWNTs, because the electrons are spread out around the circumference 
of the tube, and the probability of two electrons to be near to each other is of the 
order 1/N, where N is the number of atoms on the circumference. The interaction is 
usually screened on a length determined either by the tube length or by the distance 
to nearby gates. The Coulomb interaction causes Umklapp, backward and forward 
scattering processes among the electrons. Away from half-filling Umklapp scattering can 
be neglected. We also disregard backscattering processes, which is a valid approximation 
if the tube radius is not too small [TT]. The forward scattering processes can be fully 
included within a Tomonaga-Luttinger (TL) model for SWNTs (HI CHI HI], yielding the 
TL Hamiltonian 

H e i = H N + ^Hj , (1) 

j 
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Hj=f I da 



9j 



(2) 



where the index j = c+,s+,c— , s— counts the four excitation sectors for total charge, 
total spin and relative (with respect to the two electronic subbands) charge and relative 
spin, respectively. The first term H N = J2j EjNj /2 in Eq. (JTJ) represents the ground 
state energy. The quantities Ej are the energies for total/relative charge and spin for 
a given number of electrons N c ± and given spin N s ±. Specifically, E c+ contains the 
charging and single particle energies and is taken as a free-parameter hereafter, in order 
to include screening effects of nearby gates. For the remaining modes is Ej = Eo/4, 
with e = hvp-n/L the single particle level spacing. Here we denoted with vp the 
Fermi velocity, while L is the SWNT length. For SWNTs the term Eo5N c _ should be 
added in H N , where 5 is a small off-set caused by the discrepancy of the discretized 
energies between the two bands. The gate potential term —eN c+ Vc is also added in the 
Hamiltonian. The second term in Eq. describes the collective bosonic excitations 
of the one-dimensional interacting electron system. Hj{x) and <pj{x) play the role of 
the momentum and displacement field operators, respectively, and they are conjugate 
variables of each other, i.e., [(j)j(x), Hjt(x')] = ih5jj/5(x — x). Due to the interaction, the 
excitations propagate with renormalized velocity Vj = vp/gj, where gj is the interaction 
parameter for the j-sector. Because the Coulomb interaction describes an interaction 
among total electronic densities, only the total charge sector is affected while the other 
sectors are neutral sectors, i.e., g c+ < 1 and gj = 1 (j = s+, c— , s— ), respectively. 
Hereafter we denote g = g c+ and we call j = c+ the charge sector simply. The interaction 
parameter is roughly estimated as g = 1/yl + U/hvpn with U = 8e 2 log R sc /R. Here 
R sc is a phenomenologically introduced screening scale, long compared with the tube 
radius R. It can be estimated to be R sc = min{L, d}, where L ~ 1/xm is a typical 
SWNT length, and d is the distance to a possibly present gate electrode. The interaction 
parameter is estimated to be g ~ 0.2 for R sc ~ 1/im |10j . The screening caused by a gate 
electrode nearby the nanotube is typically more effective in screening, which pushes g 
closer to the non-interacting value, g = 1. 

We consider next the effects of the vibrations. The relevant phonon modes at low 
energies are those for which the quantized (due to the periodic boundary conditions 
in the circumference direction) transverse momentum is zero. There are three modes: 
stretching, twisting and breathing modes. In terms of the deformation potential, the 
transverse-acoustic twisting mode does not contribute at first order in the displacement 
|2()j . so that it will no longer be considered in the rest of this work. The transverse 
acoustic breathing mode has a finite frequency uj-q, of the order of (0.14/^,i?)meVA in 
the long wavelength limit. The energy scale is thus larger than that of the low energy 
electrons in the vicinity of the Fermi energy. However, as we shall discuss later, this 
mode strongly renormalizes the energy of the charge sector. In the following, we start 
by explicitly considering the effects of the coupling to the longitudinal stretching mode. 
We start from a continuum model for the vibrations |2ZI; such that the stretching mode 
corresponds to the longitudinal acoustic mode of the one-dimensional continuum. At a 
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second stage of the calculation also the effect of the breathing mode will be discussed. 
The Hamiltonian for the longitudinal phonons and the electron-phonon interaction are 
written as 



1 r L 

H L , P h = o / dx 

2 Jo 



l -P\x) + Cvi{d x u{x)? 



(3) 



L 



#ei-L, P h = c dxp(x)d x u(x), (4) 





where ( is the carbon mass per unit length, t> st is the sound velocity of the stretching 
mode, c is the electron-phonon coupling constant, and p(x) = 2d x (f) c+ (x)/\/irh+N c+ / L is 
the electron density. P(x) and u(x) are the momentum and displacement field operators, 
obeying the commutation relation [u(x), P(x')} = ih5(x — x'). 

Before we proceed, let us discuss the relation between the coupling constant c of 
our one-dimensional model and the three-dimensional nature of the electron-phonon 
coupling in SWNTs. For a SWNT the electron-phonon interaction is written as [21] 



r2nR r-L 

Hd-ph = / dy dxp(x,y)V(x,y), (5) 
Jo Jo 

where the y axis is along the circumference, p(x, y) is the density of electrons and 
V(x, y) is the deformation potential. For the stretching phonon mode V(x,y) = 
2c'd x u x p/ (B + p) has no y-dependence. Here u x is the displacement field, and B, p 
and d are the bulk modulus, shear modulus and deformation potential for a graphene 
sheet, respectively. Therefore, with p(x) = J 2nR p(x,y)dy, we have the relation 
c = 2c'p/{B + p). 

Since we want to investigate transport properties of the coupled electron-phonon 
system in a tube of finite length L, we employ open boundary conditions for the 
electrons, 11,(0) = IT,- (L) = 0, and <f)j{L) - ^(0) = PH For the phonon 

field we have similar boundary conditions, -P(O) = P(L) = 0, w(0) = u(L) = 0. The 
discreteness of the phonon and electron spectrum in a finite length SWNT is better 
visualized by expanding each field operator in Fourier series as 



^ l ^sax(k n x)-^={b\ n + b^ n ), (6) 

n>l V fc 



Ufa) 



fa ^ , — 

, —J2 sin ( k nX)ykn(bj, n - b j,n)- (7) 
\ 9j L n>l 

Here bj >n (&!•„) are annihilation (creation) operators of the j-sector's bosons, satisfying 
[bj >n , b\t n ,} = 8jj>5 nn >, and k n = ixn/L is the wave number. Similar expansions hold for 
u(x) and P(x): 



u\x) 



h l 

J2 sin(k n x)—=(al + a n ), (8) 



C^st-^ n>1 \fkr, 



p( yX ) = iJ^^J2sm(k n x)^/k~(al-a n ), (9) 
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where a n (at) are the annihilation (creation) operators of the phonon mode obeying the 
commutation relation [a n , a,] = b nr > . Using these relations we get 

H el = £ 3 nh ),J ) j,n + H Nl (10) 

#L, P h = J2B a naia n , (11) 



n>l 



#cl-L,ph = In ( a l + a n)(&tf,n + b c +,n)- (12) 



n>l 



Here £,• = e /gj is the discrete excitation energy for the j-sector, while e a = hv st 7i/L 



for the phonons. Finally, / = c^hg/Tc^v st TT/L is the electron-phonon coupling. The 
electron-phonon interaction term ()12j) reflects the fact that the phonons couple only 
to the total charge density via the deformation potential. The Hamiltonian H = 
H e \ + i?L, P h + -f^ci-L,ph is bi-linear in the boson operators characterizing the electronic 
and phonon excitations, and as such it can be exactly diagonalized with a Bogoliubov 
transformation (see Appendix). The diagonalized Hamiltonian is, 

H=J2 Epnf3i{3 n + E a nala n + £ s nb]^ c+ J^ n + H N , (13) 

n>l n>l j',n>l 

where Ep/ a comes from the branch of the charge density/phonon, defined by, 

E, 



-2 _|_ r 2 
'a + £ c+ 



± 



\ 2 2 



'C+ 




c+ • 



(14) 



Here we assumed the relation e a < e c+ , as it holds for carbon nanotubes. The summation 
j in Eq. (jTHj) is for the neutral sectors j = c— ,s+,s— . We notice that for electron- 
phonon couplings I 2 = e a £ c +/4: the energy E a vanishes. For even larger couplings E a 
becomes complex which is an unphysical situation. Thus we require that I 2 < e a e c+ / 4, 
where the equality sign defines the Wentzel-Bardeen singularity. We notice that also 
the diagonalized energies have a linear dispersion relation. Hereafter we use the index 
/i = a, P, v for the three kinds of bosonic excitations, and set E v = e . 

Let us now also include the effects of the breathing mode. A quantitative estimate 
of these effects is discussed in [20], where the authors evaluate the effective retarded 
electronic action obtained upon integrating out the phonon stretching and breathing 
modes. It turns out that retardation effects can be disregarded for the breathing but 
not for the stretching mode. Thus the breathing mode yields a renormalization of the 
velocity associated to the charged plasmon mode as v* + = a~BV c+ , where 



a * = f~J[9 2 i ( 15 ) 

and Rb — 2.4±0.9A j2U]- With v f = 8 x 10 5 m/sec, v st = 2 x 10 4 m/sec, and interaction 
parameter g = 0.2 the charge velocity v* + becomes comparable to the sound velocity t> st 
of the stretching mode if — 5 x 10~ 3 . This cannot occur, not even for the smallest 
SWNT radii. However, the renormalized charge velocity might become of the order 
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of the sound velocity v st for realistic nanotube radii for strongly screened interaction 
(yielding values of g close to one). The breathing mode also yields a renormalization 
of the electron-electron interaction parameter as g* = vp/v* + = g/a&- This results 
in the increase of the electron-phonon coupling in Eq. (|12|) with (J*) 2 = I 2 /a. In 
the following, we shall then consider breathing-mode renormalized parameters, with a 
parameter a-Q = 0.993 (g = 0.2), 0.804 (g = 1) as for a realistic (10, 10) tube. 

We notice that also Umklapp scattering processes should reduce the charge velocity 
|29j . This effect would be dominant near half-filling, where a strong coupling with 
the stretching mode is expected |28] . Moreover, a suppression of the electron-phonon 
coupling is expected at half-filling |3TT| 13 lj. This situation has not been considered in 
our calculation, since SWNTs are usually away from half-filling [llj . 



3. Phonon assisted tunneling 

With the help of Eq. (fTB^) we can now evaluate the density of states p(e). In particular, 
we shall focus on the density of states at a nanotube end p e nd, yielding direct information 
on the differential conductance. The density of states is calculated by expressing the 
electron field operators in terms of the boson operators [01 EH ITTj . After lengthy 
but standard calculations, see Appendix, we arrive at the following finite temperature 
density of states at the end of the nanotube: 

Pcnd(^) = J! X] X Ca,lC{3 tm C Ut n 
s=± A,cr l,m,n 

x 5 (hu - s{E s Xa + lE a + mEp + nE v j) , (16) 

where s is the index for the positive and negative energy region, A = ±, a — ± 
are the band and spin indices, the summations l,m,n run on all integers. E Xa = 
J^jCjEjNj + XeoS — eVc and = E c+ + 3eo/4 — E Xa , where Cj = l,X,a,Xa for 
j = c+, c— , s+, s— , respectively. When the temperature is much lower than the level 
spacings, k^T E^, the coefficients C^ p are written as 

C^ p =C^ p + 7 At e E ^/ k ^ T (C^p+i + — 2C^ p °) , (17) 

c -° = T -M^ c ^° e(p) - (18) 

Here T is the Gamma function and G is the step function. The coefficient 7 M corresponds 
to the end-tunneling exponents in the bosonic correlation functions, with 

1 E a 9 1 Ea o 

7a = ~, sin cp, 7^ = cos ip, (19) 

4^ e c+ Ag e c+ 

and j v = 3/4. The phase cp is defined by <p = arctan(— AI ^/e a e c+ / (e 2 a — e1 + ))/2. 
Thus, from (jl6j). one sees that at low temperatures a signature of the finite length 
of the SWNT appears in the form of discrete peaks in the density of states. More 
precisely, three peak families with distinct periods are identified, cf. Fig. 1, reflecting 
(i) the discrete excitations of the (renormalized) charge density with energy spacing 
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Ep, (ii) the excitations of the neutral sectors with spacing E u , and (hi) the phonon- 
assisted-like excitations with spacing E a . At zero temperature the l-th phonon peak 
height with respect to that of the corresponding electron excitation peak is given by 
CaJ°/Cafi° = 1a(la + 1) • • • (7a + I — 1)/U. For example, the ratio for / = 1 is simply 
7 Q . For large 7^ it holds the asymptotic relation 

C l? - |^ T o°©b) , (20) 

which resembles the relation among satellite peaks due to vibron-assisted tunneling as 
in [32] ■ For large p, on the other hand, each peak series shows power-law in its height, 
because the relation C^=° oc p" 1 ^ 1 holds if one uses the asymptotic expansion for the 
Gamma function. 




frco [meV] ha [meV] 



Figure 1. End-tunneling density of states of a (10,10) SWNT vs. frequency for 
strongly interacting electrons (g = 0.2, left figure) and non-interacting electrons (g = 1, 
right figure). Only the positive energy region is shown. We use the parameters 
E+ + = E+_ = 2mcV, E+ + = 2.99meV, and E+_ = 1.34meV for both figures. This 
parameter set corresponds to the condition iV c _ = — 1, N s+ = 1, iV s _ = — 1, and finite 
band off-set 5 = 0.2, in which there is an unpaired electron in the A = — band so that 
a spin-1/2 appears in the SWNT. The no peak region in < hui < 1.34 is the Coulomb 
blockabe region. The large three peaks at 1.34, 2 and 2.99mcV reflect the different 
addition energy to each level in the SWNT. For each of these addition energy peaks, 
three families of peaks are clearly identified. They are associated to the plasmon, 
phonon and "neutral" bosonic excitations of the suspended SWNT Ep, E a and E„, 
respectively. Because the relation E_, = E^_ + E v holds for the present electron 
configuration, some of the excitation peaks stemming from fnu = E^_ are degenerate 
with some of those stemming from hid = E^ + . The remaining parameters for the 
left figure are E a = 0.0410meV, Ep = 8.23meV, E v = 1.66meV, -y a = 1.52 x 10~ 9 , 
jp = 1.24, 7„ = 3/4, T = 0. For the right figure, E a = 0.0327meV, Ep = 1.33meV, 
E v = 1.66meV, 7 Q = 1.76 x 10~ 6 , 7,3 = 0.201, -y„ = 3/4, T = 0. 



Here we estimate the amplitudes of the phonon peaks and the excitations' spacings. 
First, we use the following quantities (TQIEII: v f = 8 x 10 5 m/sec, v st = 1.99 x 10 4 m/sec, 
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c = 20eV, R = 6.79A (for a (10,10) armchair SWNT), ( = 2itRM, (M = 3.80 x 
10~ 7 kg/m 2 ), and L — 1.0 x 10~ 6 m. These values give (for g = 0.2 then a B = 0.993) 
the energies E a = 0.0410meV, Ep = 8.23meV, E v = 1.66meV, 7 Q = 1.52 x 10~ 9 , and 
■jp = 1.24, The estimated peak ratio 7^ is quite small. However, we note that the value 
is sensitive to the parameters. We may have g ~ 1 for the screening caused by the gate 
electrode. For g = 1 (a B = 0.804), E a = 0.0327meV, Ep = 1.33meV, E v = 1.66meV, 
7 a = 1.76 x 10~ 6 , and 7^3 = 0.201. The sensitivity of 7 a will be discussed in the next 
paragraph. We also note that there is an uncertainty in the precise value of the coupling 
constant c. Estimates in [27] yield values c ~ 14 — 21eV. 

In Fig. the ratio of the first phonon peak height 7 a as a function of the charge 
plasmon energy spacing and of the electron-phonon coupling is plotted. For very strong 




Figure 2. Relative phonon peak height ~f a as a function of the ratio between charge- 
plasmon and longitudinal phonons level spacing e c+ /e a , and of the electron-phonon 
coupling I/e a . The vertical axis is in logarithmic scale. The back side region without 
mesh is the WB singularity region, where j a has an imaginary value. 

couplings, corresponding to the condition I 2 > e a e c+ /4:, the system becomes unstable, 
as the Wenztel-Bardeen (WB) singularity is approached. In general the ratio is very 
sensitive to the parameters. The peak height becomes larger for small charge velocities, 
especially close to the WB singularity. For the small coupling case, the system survives 
the collapse for a wide range of parameters while the charge velocity is increased, until 
finally the approach of the singularity is expected. 

4. Conclusion 

From our semi-quantitative discussion, it follows that the phonon peaks could be 
observed in experiments. In recent experiments on suspended metallic carbon 
nanotubes, many excitation peaks have been observed in the differential conductance. 
The peaks appear outside of the Coulomb diamonds, and are almost equidistantly 
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spaced. The period of the peaks is about 0.53meV for a L = 140nm length nanotube, 
0.69meV for a L = 420nm, and 0.14meV for a L = 1200nm, respectively. The theoretical 
estimate of the period for the phonon assisted peaks is 0.77 to 0.09meV, while that for 
the charge density is 47.5 to 5.5meV (if g ~ 0.2), and 9.5 to l.lmeV for the neutral 
sectors. The observed peaks might be the phonon-assisted peaks discussed in this paper. 
The relative peak height 7 Q reflects the strength of the electron-phonon coupling. From 
the experimental data in Ref. [33] a rough estimate 7 a ~ 0.3 to 0.8 can be extracted. 
These values are much larger than those we estimated in the previous section for strong 
Coulomb interaction (g ~ 0.2). Moreover, the experiments seem to well agree with the 
asymptotic formula (|2T)|) . This could be an indication of a strong renormalization of the 
electronic parameters expected when the system is close to the WB singularity. Thus, 
while a qualitative agreement with the experiments is found, an understanding of the 
lack of a quantitative agreement clearly calls for further theoretical investigations. We 
mention that a relation similar to ()20|) has also been predicted for electronic degrees of 
freedom dressed by a Frank-Condon factor originating from a generic vibron-electron 
coupling. However, also the theoretical expectation for the values of the coupling in the 
model in seems to be much smaller than the one extracted from experiments |33j . 

For temperatures larger than the level spacing finite size effects are washed out, 
and we can adapt our calculation to the case of a SWNT of infinite length. The density 
of states for this case becomes: 



where g' is the interaction parameter renormalized by the electron-phonon interaction, 
(g 1 ) -1 = 4(7 a + 7/3). As expected, the charge, neutral and phonon peaks are washed out 
and the density of states shows a simple power-law behavior with exponent modified 
by the electron-phonon interaction. The renormalized parameter g' is always larger 
than g, showing that the phonons contribute to an effective attractive electron-electron 
interaction. 

When a gate voltage is applied to the nanotube, the possibility to investigate 
Coulomb blockade phenomena and Kondo-correlation effects arises. In the presence 
of a gate potential the coupling between the electrons and the transverse bending 
mode should be included. One can control the coupling strength by changing the gate 
potential. Because the bending phonon mode couples to many of the charge density 
modes, and it has less energy for the vibration, an instability might be easily induced. 
The effects of the gate potential will be discussed elsewhere. 

In conclusion, we have discussed transport properties of a suspended metallic 
SWNT. The low energy excitations caused by the stretching phonons are expected to 
appear as small periodic peaks beside the peaks for the charged and neutral electronic- 
like modes in the differential conductance. Each series of peaks exhibit a different 
power-law. When the system is close to the WB singularity, the phonon peaks are large 
enough in height to be observed. The observation of the phonon peaks would thus 
indicate that the singularity is approached. 




(21) 
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Appendix: Density of states of suspended SWNTs 

Here we show the detailed calculation of the density of states at the end of a finite 
length nanotube. Let us consider a high enough tunneling barrier between electrodes 
and the nanotube. At low temperature, the differential conductance is proportional to 
the density of states at the end of the nanotube 

dI/dV hiBS oc p C nd(^ = eVbias) = -lmG cnd (hu = eV^A , (22) 

where G en <i{Tiu)) = G(x,x';hu>)\ x=x >^ is the Fourier transform of the retarded Green's 
function G(x, x'; t) = ±9(t)(V(x, t)&(x', 0) + &(x', 0)*(x, t)), and f is a length which 
is of the order of the inverse of the Fermi wave length. The electron field operator is 
written in terms of boson operators using the bosonization identity 



* v (x,t) = —. = u v e iK i x e-^ Cj6dAx ' t \ (23) 
V2vre 

and *&(x, t) = J2r]^r)(x,t). Here rj = (r, A, a) is the parameter set of the left and right 
going wave (r = — ,+), the band index of the nanotube (A = ±), and the spin index 
(a = ±). e — > is the cut-off length, U v is the Klein factor, k v = r\k F + r(N v + XS)tt/L 
and N v is the number of electrons in the ^-branch. k F = 2ir /3a is the Fermi wave length, 
a is the honeycomb lattice constant, and we consider the case hp 3> iN^n/L in this 
paper to justify our low-energy theory. The field 6j jT is written as, 

hr(x, t) = (-r4>j{x, t) + J X dx%{x')\ . (24) 

We note that the left and right going waves are not independent anymore in the case of 
open boundary condition; the relation f r ^ ](r (i) = — A CT (— x) holds. 

Let us consider the correlation function in the Green's function t)W{x', 0)). 

This can be separated into the boson and the fermion parts; 

(¥(z,t)*V.0)>= £^W)^(°)> 

x ^ e i ^v{ x - x ')Y[.^Q- i \ c j e jA x ^)Q i \ c 3 e 3A x 'fi)^ 

_ Q iK v( x + x ')]J^ e - i ^ c 3 9 j.r( x ,t) e i^C]9j : r(- x \0)^ ^ 

Since the Klein factor only changes the number of electrons in ^-branch, the fermion 
part of the correlation function can be calculated as, <C/„(t)£/}(0)) = exp(-i Nri+ \ Nr > t). 
Here Ejy + \ — En v = E^ a being addition energies of the 7y-branch. (Because of the 
open boundary condition, we have no r dependence in the addition energies). On the 
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other hand, using the relation for bosonic operators, (e tA e tB ) = e a^ 2 ) \( b2 )+( ab ) ) the 
bosonic part of the correlation function can be written as 



(e- 



,{-i(e?(x,t)>-i(« r («',o)>+(fl 3> (x I t)e 3 -, r (x',o)>}/4 



(26) 



which doesn't depend on A, a. Therefore the calculation can be reduced to that for 
the bosonic correlation function, (9j jr (x, t)9j ir (x', 0)}. Now we discuss the correlation 
function of the charge sector. Similar calculations have been done for the neutral sectors. 
Substituting the Eqs. ©, (J7J) into (JUJ), an d using the Bogoliubov transformation 



3n 

at 

V Ct n j 



ARC 

r ^n- LJ n K - y n 



J c+,n 
V a n j 



(27) 



where 



A 



'E n 
1 



-ijEpn 



'Ean 










i^jEpn 

^E~H 









/ E a n 



(28) 



( COS if 




V sunp 



e c+ n 
i 







vfc y 

— sin ip \ 

cos ip smip 

- sin if cos y> 

cosy? / 





1 i 



(29) 



2\/ e c+ n 
1 

2 V /£ C+ 71 













(30) 
(31) 



at the end of nanotube x 



x 



£, the correlation function satisfies the relation 



(9 c+tr (x,t)9 c+tr (x',0)) 



1 7T 



9 3 L n>l K ri { £c+ 



+ 



sin 2 v9(a n (t)«t (0) + a£(f)a: n (0)) 



t, 



(32) 



■-c+ 



which doesn't depend on r. The lowest order of the function with respect to £ is 
therefore we calculate the bosonic correlation function at x = x' = hereafter. Using 
the relation 

((3 n (t)(3l(0)+f3i(t)(3 n (0)) = cos(^nt) cosh -isin(^nt), (33) 
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and multiplying by e~ hvkn ^ Ecut each term of the summation over n in Eq. (jH2| to 
avoid divergences in the calculation, where E cnt = Tive^ 1 ^-^ oo) is a cut-off energy 
and v = E a /pL/7rh, then the exponent of the Green's function is written as 

(9l r (0,t)) + (9 j;r (0,t)9 j>r (0,0)) 

l E f l _ e -Ep/E cute -iE t/h 

= -^r os n log — — 
. Vlflg ^Ma+0)-cos(ft) 

cosn^ cut + fcBr J i 
+ {/5 — > a, cosy? — > sin<£>}, (34) 

where we use the identities 



coshx =l + 2^e- 2nx , (35) 

n=l 

oo „—nx 

E = -log(l-e-). (36) 

Then the boson part of the correlation function is calculated as 

(p-faAW^iAW)) = Da (t)Dp(t), (37) 
where, up to orders (exp{— E^/hsT}), 



Dp(t) 



c 



cut 



"7/3 



x|| , ^(g; + ^)-cos(fO 
cosh ( rz: + i^t) - 1 



E A^"^, (39) 



p=— oo 



where 



^ = ^° + 7/3e"^ /fcBT (l^ + - 2D™) , (40) 

and (3 — > a for D a (t). After a similar calculation for the neutral sectors and for 
(x' , 0)^ (x , t)) , we finally get the Green's function 

G en d{hu>) = — - — - — — E E E D at iDp m D vn 

s=± A,cr l,m,n 

x (42) 

hto - s(E s Xa + lE a + mEp + nE v ) + i0 + ' v ' 

After re-definition of the coefficients, e.g., C^ p = -D M , P (2 sin 2 (/cF£)/ 7re ) 1//3 ; we § e t the 
density of states ([TBI) . 
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